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CUBIC HAMILTONIANS 


P.L. ROBINSON 


Abstract. We determine a precise necessary and sufficient condition for completeness of 
the Hamiltonian vector field associated to a homogeneous cubic polynomial on a symplectic 
plane. 


0. Introduction 

The flow of the Hamiltonian vector field generated by a smooth function on a symplectic 
manifold is a familiar object of study. Let the symplectic manifold be simply a symplectic vector 
space: the Hamiltonian flow generated by a homogeneous linear function is a one-parameter 
group of translations; the Hamiltonian flow generated by a homogeneous quadratic function is 
a one-parameter group of linear symplectic transformations. In each of these two cases, the 
Hamiltonian flow is complete: each maximal integral curve of the Hamiltonian vector field is 
defined for all time. The case of cubic Hamiltonian functions is different: for some cubics the 
flow is complete whereas for others it is incomplete. 

Our primary objective in this paper is to establish a simple necessary and sufficient condition 
for the cubic "0 on a symplectic plane (Z, H) to generate a complete Hamiltonian flow. In Section 
I we associate with ip a suitably symmetric linear map from Z to the symplectic Lie algebra 
sp(Z, H); following this map with the determinant yields a quadratic map A : Z —> K. In 
Section 2 we analyze an arbitrary integral curve z : I Z oi the Hamiltonian vector field 

OO 

defined by ip] we find that the second time-derivative 2: equals 2 Fz, where the scalar function 

°° 

F:=Aoz:J—>K satisfies the equation F = QF"‘ familiar from the theory of elliptic functions. 
In Section 3 we achieve our primary objective, proving that the Hamiltonian vector field 
is complete if and only if the determinant A is identically zero; beyond this, we comment on 
the nonconstant integral curves of in the complete case and the incomplete case. Finally, 
we assemble several remarks on issues arising from the main body of the paper: in particular, 
we remark that A is identically zero if and only if 1/; is a monomial; these remarks we plan to 
develop more fully in subsequent papers. 

In a subsequent paper we also plan to present a similar treatment of quartic Hamiltonian 
functions; for now, we merely note one difference between the cubic case and the quartic case. 
In the cubic case, the scalar function F satisfies the differential equation F = whose elliptic 
solutions are always Weierstrass Pe functions associated to triangular lattices, with 52 zero; in 
the quartic case, the corresponding scalar functions include Weierstrass functions associated to 
rectangular lattices, with 52 nonzero. 


1. Symplectic Algebra 

Let (Z, H) be a real symplectic vector space: thus, Z is a vector space and H : Z x Z —>■ R a 
nonsingular alternating bilinear form. Though it is not necessary for some of what we shall say, 
we suppose throughout that Z is two-dimensional, so that (Z, H) is a symplectic plane. The 
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symplectic algebra sp(Z, fl) is the (commutator bracket) Lie algebra comprising all linear maps 
C : Z ^ Z such that for all x,y € Z 

n{Cx, y) + f2(x, Cy) = 0. 

As a vector space, sp(Z, fl) is canonically isomorphic to the space of all symmetric bilinear 
forms on Z: to C G sp(Z, fl) there corresponds the symmetric bilinear form 

Z X Z ^ M. : {x,y) 1 -^ n(x, Cy). 

Now, let i/j : Z —>■ R be a homogeneous cubic polynomial. To ip we associate the (fully) 
symmetric trilinear function ZxZxZ^M. with value a.t {x,y, z) € Z x Z x Z given by 

'^{x,y, z) = ip{x + y + z) - {ip{y + z) + ip{z + x) + ip{x + y)} + ip{x) + ip{y) + ip{z). 

When z € Z is fixed, 'I'(x,t/, z) is symmetric bilinear in {x,y) € Z x Z; it follows that there 
exists a unique G sp(Z, n) such that for all x,y G Z 

'^ix,y,z) = 2il{x,rzy). 

Full symmetry of guarantees that the resulting linear map 

F’^ = r : Z sp(Z, n) 

is symmetric in the sense that for all x, y G Z 


Note that if z G Z then 


F^y - T'yX. 


2n{z, r^z) = ^'(z, z, z) = {27 -(3x8) + 3}ip{z) = 6ip{z) 


or 

= ^f^(^,r;2Z). 

Differentiation of this formula for ip yields the result that if i;, z G Z then 
V'z(v) = i{D(v, F^z) + 0(z, F„z) + D(z, F^n)} 

whence by symmetry of F : Z —>■ sp(Z, D) it follows that 

ip'^{v) = n{v,T^z). 

Of course, as ■;/' is a cubic, the first derivative ip'^ is quadratic in z G Z. As a bilinear form, the 
second derivative ip” at z G Z furnishes another means of introducing ik and F: indeed, if also 
x,y G Z then 

i^ziy, x) = ^'(x, y, z) = 2n{x, T^y). 

This equation represents ip” by 2 F 2 relative to the symplectic form D; consequently, the classical 
Hessian of ip is Det(2F2). 

According to the Cayley-Hamilton theorem, if z G Z then 

r,r, - (Tr F,)r, + (Det r,)i = o 
whence the fact that F^ G sp(Z, D) is traceless implies that 

F.F, = -(Det r,)I. 

We define the scalar function A'^ = A : Z ^ R by requiring that for each z G Z 

A(z) = -(Det F,) 

so that 

F.F, = A(z)/. 
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Theorem 1. If z Z then A(rzz) = A(z)^. 

Proof. If 2 ; = 0 then both sides of the alleged equation plainly vanish. If z ^ 0 then apply the 
special case = f'zf'zZ of symmetry repeatedly: a first application gives 

A{r^z)z = Tr^zf'r^zZ = Tr^zI'zf'zZ = Tt^zA{z)z 
and a second application gives 

A{z)Ty-^zZ = A{z)TJ'zZ = A{z)A{z)z = A{zf'z 
whence the alleged equation follows by cancellation. □ 

2. Cubic Hamiltonians 

We shall now view (Z, H) as a symplectic manifold in the natural way. Thus, the vector 
space Z is naturally a smooth manifold; if z & Z then there is a natural isomorphism from the 
vector space Z to the tangent space TzZ sending v € Z to the directional derivative operator 
v\z € TzZ given by the rule that whenever / : Z ^ K. is a smooth map, 

v\z{f) = f'z{v) = ^f{z + tv)\t=o. 

Also, H serves double duty as a nonsingular alternating bilinear form on the vector space Z 
and as a nonsingular closed two-form on the smooth manifold Z; explicitly, if x,y, z € Z then 
the value of the two-form at 2 : is given by 

flz{x\z,y\z) = ^{x,y). 

When / : Z ^ K. is a smooth (Hamiltonian) function, the corresponding Hamiltonian vector 
field € VecZ on Z is defined by the requirement 

n = -df 

where j signifies contraction as usual. An integral curve of the vector field is a smooth map 
2 : : / —>■ Z (on some open interval / 9 0) satisfying the Hamilton equations: for each t € I the 
tangent vector to 2 : at t equals the value of at Zt, thus 

° ef 

We shall focus on the case of a homogeneous cubic i/; : Z —> K. as Hamiltonian function. The 
value of at 2 : € Z is a vector made tangent at z: say 

^t=X^{z)\z 

with : Z —> Z a smooth vector-valued function. Now, let w, 2 : € Z: on the one hand, 
(5 ^lH),(uU) = nzift,v\z) = nz{x^iz)\z,v\z) = nix^izfv); 
on the other hand, 

-difzivlz) = -ip'zix) = -fl{v,rzz) = ft{TzZ,v). 

As the symplectic form H is nonsingular, it follows that 

x'^{z) = TzZ. 

Accordingly, the Hamilton equation for 2 : : / Z reads 

Let 2 : : / —>■ Z be a solution of this Hamilton equation. Take a further derivative: as L is 
symmetric, 

2 : = ro2: -i- Lz 2 ; = 2rz 2 : = 2rzrz2: 

by a further application of the Hamilton equation. Recall that if w € Z then = A(w)I 

and write 

F := Ao z : I 
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It then follows that z \ I ^ Z satisfies the second-order equation 

OO 

z = 2Fz. 

Note here that A is defined on the whole space Z while F is defined only along the integral 
curve z. 


Theorem 2. The scalar function F satisfies the second-order equation 


F = 6F 


2 


Proof. From the definition 


FI= 


we deduce by repeated differentiation that 


FI = ror^-hr^ro 

2 2 

and 

OO 

F I = FooF^ + 2 roro -h F^Foo. 
2 2 2 2 

Here, the first and last terms on the right both equal 2 FrzF 2 
while FoFo equals F^I on account of z = F^z and Theorem [T] 


Q OO 

2F I on account oi z = 2Fz 

□ 


We may at once deduce a first-order integral of this second-order equation: multiply through 

O 

by 2 F to obtain 

O OO o ° 

2 FF= 12F^ F 

from which there follows 

(Ff =4^3-33 

for some real constant 53 . This notation is deliberately chosen to accord with the theory of 
elliptic functions. In fact, the solutions to this first-order differential equation are as follows: 

• if g 3 is nonzero then F{t) = p{t — a) for some real a where p is the Weierstrass Pe function 
associated to a triangular lattice (the so-called equianharmonic case); 

• if g 3 is zero then either F{t) = {t — a)“^ for some real a or F is identically zero. 

Note that when F is a (shifted) Weierstrass Pe function, z = 2Fz is a (vectorial) Lame 
equation and may be solved accordingly; for example, see page 285 of [Forsyth]. 


3. Completeness Characterized 

We continue to let F : Z ^ sp(Z, H) be the symmetric linear map corresponding to the 
homogeneous cubic ^ : Z —>• R on the symplectic plane {Z, H); we also continue to let z : / —>■ Z 

be an integral curve of the associated Hamiltonian vector field We shall suppose that the 

. ® 

curve z has initial point zq and hence initial velocity zq = Fz^Zq. Our aim in this section is 
to decide precisely when such an integral curve may be defined for all time; that is, precisely 
when the maximal domain of definition J is R. itself. 

The critical case is decided immediately. Let (equivalently, dip) vanish at zq; thus, z has 

O 

initial velocity Zg = Fz^zq = 0. In this critical case, the solution z : I ^ Z is plainly given by 
Zt = Zq for all t G / and the maximal / is indeed M. In this connexion, note further that if an 
integral curve z : I ^ Z vanishes at any point then so does its velocity vector and hence z itself 
is identically zero. 

Now let the integral curve z : I ^ Z he other than critical: thus, FzpZo = Zg yf 0 and of 
course zg yf 0. We distinguish two cases. 
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For the first case, suppose there exists some s G I such that 0 ^ F{s) = A(zs) and therefore 

oo o 

F {s) =F {sf > 0. The comments after Theorem [5] show that F has a double pole at some real 
a; thus = F(t)I is unbounded as t ^ a and so zt itself is unbounded as t —>■ a. In this 

case, the maximal domain of 2 ; omits a and thereby falls short of R. 

For the second case, suppose that F{t) =0 whenever t G I. Note that the linear map F^^ 
kills F^gZo (because F^^F^^ = F{0)I = 0) but does not kill zq (because F^^zq = ^0 0); thus 

O 

Zq and Zq constitute a basis for the plane Z and so 

{s(zo t Zq'J s^t ^ IR.} = i^Z \ M 2^0) U { 0 }. 

00 o 

The supposition F = 0 implies that 2 : = 2Fz = 0 so that zt = zq +1 Zq for alH € /; essentially 

O 

as in the critical case, the maximal I is therefore R. Now A vanishes on zq + ^ whenever 
t G M. (as F is identically zero) and hence vanishes on s(zo + t zq) whenever s,t G R (as A 
is homogeneous); the continuous function A now vanishes on the dense set (Z \ R zq) U {0} 
and therefore vanishes on the whole of Z. This proves that if A vanishes on the image of some 
non-critical integral curve then A vanishes identically. 

We may now marshal these facts towards the following result. 

Theorem 3. Let tp : Z ^ M. be a homogeneous cubic and A'^ the associated determinant. 

• If A’^ = 0 then is complete; each non-constant integral curve is an affine line. 

• If A'^ ^ 0 then if is incomplete; only the constant integral curves are defined for all time. 

OO 

Proof. If A = 0 then each maximal integral curve z has F = 0 so that 2 : = 2Fz = 0 and z on 
R is affine, as we have seen. If A ^ 0 and the integral curve z is not critical, then F ^ 0 so 
that z experiences finite-time blow-up, as we have seen. □ 

Looking ahead to the next section, we remark that A’^ is identically zero if and only if ip is 
monomial in the sense that there exists w G Z such that for all z G ^ 

4. Remarks 

In this closing section, we record a number of miscellaneous remarks that stem from the 
body of this paper. 

COORDINATE EXPRESSIONS 

Though our whole approach has been intentionally coordinate-free, it is also of interest to 
see the development in terms of linear symplectic coordinates, not least because this may offer 
glimpses of a fresh perspective on classical invariant theory. 

To this end, let u,v G Z satisfy Ll{u,v) = 1 and so constitute a symplectic basis for {Z,Ll). 
Decompose z G Z as 

z = pu qv 

with 

p = p{z) =n{z,v), q = q{z) = n{u,z). 

Write 

a = n{u, F„u), b = Ll{u, F^m), 
c = Ll{v,Tuv), d = Ll{v,Tyv). 

With these conventions, the cubic 

^’( 2 ) = ■^f^(2:,r22:) 
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has coordinate form 


-tpiz) = + ibp^q + icpq^ + dq^} 

o 


and the (vector) Hamilton equation 
becomes the familiar scalar pair 

The associated determinant 


z = T,z 


o dib o dtb 

P = -9 = —. 

dq' dp 


A(z) = -(Det r,) 


assumes the form 

A(z) = {b^ — ac)p^ + {be — ad)pq + (c^ — bd)q^ 
and is the Hessian of ip (up to scale). We are not the hrst to observe that the discriminant 

{be — ad)'^ — 4(6^ — ac)(c^ — bd) 
of this quadratic is precisely the discriminant 

a^d^ — 3b^c^ + 4ac^ + 46^d — 6abcd 


of the cubic 

ap^ + 3bp^q + 3cpq^ + dq^; 
for example, see page 60 of [Salmon]. 

Of course, a purely coordinate-based approach is possible. Let us indicate partial derivatives 
more succinctly by means of subscripts. With the cubic 

ip{z) = -f Sbp'^q + 3cpq^ + dq^} 

o 

as above, direct computation reveals that ippqipq — ippipqq is divisible by p and ipqpipp — ipqippp is 
divisible by q; in each case, the quotient is precisely 2{{b'^ — ac)p‘^ + {be — ad)pq+ {c'^ —bd)q'^} and 
we recover (twice) the determinant A in coordinate form. In fact, when the Hamilton equations 

P=-i^q. d=1pq 

are differentiated by time once more, they yield precisely 

oo oo 

P ~ '^pq'^q '4^pd^qqi ^ ~ '^qp'4^p '4^q'^pp 

oo 

and we recover the scalar components of z = 2Fz. 

CANONICAL FORMS 


The simplest type of homogeneous cubic is a monomial: for w G Z define ip'^ : Z —> R by 
requiring that for all z € Z 

ip-{z) = ^n{w,z)\ 

For this cubic, the corresponding symmetric linear map T™ : Z -A sp(Z, H) is given by 

= fi{z,w)n{w,v)w 


whenever z,v G Z, and the associated determinant A™ is identically zero. 
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Conversely, let the cubic with corresponding symmetric linear map T be such that the 
associated determinant A is identically zero. We claim that i/; = '0™ for a unique w G Z; to 
justify this claim, we may of course assume that F is not itself identically zero. Note that if 
z G Z then F^F^ = 0 so that Ran F^ C Ker F^ with equality precisely when F^ ^ 0. Note also 
that if x,y G Z then 

Fa;Fj^ + FyFj; = {A(x + y)- A{x) - A{y)}I = 0. 

When x,y, z G Z let us write 

liX,y,Z) = T^VyZ. 

Observe that this expression is now antisymmetric in its first pair of variables and was already 
symmetric in its last pair; thus 

y, z) = jix, z, y) = - 7 ( 2 ;, x, y) = --f{z, y, x) = 'y{y, z, x) = 7 ( 1 /, x, z) = - 7 ( 0 ;, y, z) 
and so 7 vanishes identically. This proves that if x^y G Z then 

Ran Fy C Ker F^, 

and choosing any z G Z with F^ 7 ^ 0 then gives 

Ran F^ C Uyg^Ran Fy C na,g 2 Ker F^, C Ker F^ 
with equality of the end terms and hence equality throughout, whence 

Uyg^H-^Ti Fy 03 ,g^Ker F^. 

is a distinguished line in the plane Z. Let w G Z he & basis vector for this line. If z G Z then 
Tz = Xzi-)w for some linearly z-dependent Az in the dual Z*: as F^ kills w so does Az and 
therefore Az = ■) for some /iz G K also linear in z; this shows that 

F^ = yz^iw, ■)w 

for some y G Z*. Symmetry of F forces y to kill w so that y = un(-, w) for some u € R.. In the 
resulting formula 

Fz = h'n(z,w)n{w, ■)w 

the cube root of the scalar v may be absorbed into w; this renders w unique and we conclude 
that F = F“ as claimed. 

Thus, the assignment w i-A F“ is a (cubic!) bijection from Z to the set of all symmetric 
linear maps Z sp(Z, fl) for which the associated determinant A is identically zero. 

The same conclusion may be reached efficiently (though prosaically) using coordinates. From 
the identical vanishing of A in the form 

(&^ — ac)p^ + {be — ad)pq + (c^ — bd)q^ = 0 

we deduce (by setting y = 0,p = 0, and pq^hiu turn) that 6^ = ac, = bd, and ac = bd. Let 

A be the cube root of a and y the cube root of d: then 

(X'^y)^ = a^d = a ■ ad = a ■ be = b ■ ac = b ■ b^ = b^ 

so that y = b and Xy^ = c likewise; it follows that the cubic is a monomial, namely 

ap^ + Sbp^q + Hcpq^ + dq^ = {Xp + yq)^ ■ 

When the determinant A is not identically zero, there are three possibilities: 

• A(z) = 0 for 2 ; on a line-pair through 0 and A takes values of each sign elsewhere; 

• A( 2 :) = 0 for z on a line through 0 and A is positive elsewhere; 

• A(0) = 0 and A is positive elsewhere; 

and canonical forms may be developed for each of these. In connexion with these possibilities, 
we remark (from Theorem [T]) that if A takes negative values then it also takes positive values. 
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EVALUATION OF ga 

Let tp : Z ^ M. he a, homogeneous cubic and let the Hamiltonian vector field have 

OO 

0 : / —>■ Z as an integral curve. As we have seen, z = 2Fz where the scalar function F : / —s- R 

O 

satisfies (f)^ = 4F^ — 53 for some constant 53 that depends on the integral curve z. 

Let the initial point Zg be such that ^^(zo) = 0; as the Hamiltonian ip is constant along the 
integral curve, it follows that ip{zt) = 0 for all t G I. If Zq itself is zero, then of course F = 0 
and gs = 0. Now assume that Zq is nonzero, so that Zt is nonzero for all t € I. For each t G I 
we have 0 = 3ip{zt) = H(zt, Zt) whence (as Z is a plane) Zt is parallel to Zt; say z= X z for some 
scalar function A : / —?■ R. On the one hand, 

Fz = = Fz z = FzAz = AFzZ = A z = A^z; 


on the other hand, 
Thus 


2Fz — z — A.^FAz — X z A^ 


A z + Fz. 


X = F = X 


2 


and so 

F = (a2)° = 2A a = 2AF = 2A^ 

It follows that in this case. 


33 = 4F3 - {pf = 4{X^ f - {2X^f = 0. 

In short, an initial point Zq with ip{zo) = 0 spawns an integral curve for which 33 = 0. 

Let us offer some sample computations in coordinates, li ip = — q^) then P = —ipq = 

000 OO o 

and q = ipp = so that P = 2{pq)p and q = 2{pq)q; thus F = pq so F = Fqipp — Fpipq = p'^A-q'^ 
and 33 = 4F^ — (f)^ = Ap^q^ — {p^ + = —{p^ — q^Y ff 3 = ^ 0. Similarly, if 

O 

Ip = p'^q + pq^ then F = p^ pq -\- q^ and F = (q — p)(2.p + q){p + 2q); after considerable 

O 

simplification, 33 = 4F^ — (f)^ yields 33 = 27ip'^ > 0. 

Finally, we remark (without proof - but see page 100 of [Salmon]) that classical invariant 
theory reappears in general: if 

6 = + 4ac^ + Ah^d — habcd 


denotes the discriminant of the cubic Sip then 

33 = -9S ip'^ 


SO 

{pf =AF^ + 95ip‘^. 
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